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Based on the non- autonomous quantum master equation, we investigate the dissipative and de- 
coherence properties of the two-level atom system interacting with the environment of thermal 
quantum radiation fields. For this system, by a novel algebraic dynamic method, the dynamical 
symmetry of the system is found, the quantum master equation is converted into a Schrodinger-like 
equation and the non-Hermitian rate (quantum Liouville) operator of the master equation is ex- 
pressed as a linear function of the dynamical u(2) generators. Furthermore, the integrability of the 
non-autonomous master equation has been proved for the first time. Based on the time-dependent 
analytical solutions, the asymptotic behavior of the solution has been examined and the approach 
to the equilibrium state has been proved. Finally, we have studied the decoherence property of the 
multiple two-level atom system coupled to the thermal radiation fields, which are related to the 
^"V quantum register. 
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I. INTRODUCTION 

The behavior of the quantum dissipative systems interacting with the background thermal radiation field, is one 
of the central subjects in quantum statistical physics. Extensive interests in these systems arise from many fields 
£> \ of physics, e.g., condensed matter physics Q], quantum optics 0, 03, H, quantum measurement 0, quantum 
■ computation 0, and so on. The dissipation and decoherence are generated due to the interaction between the 
system and the thermal bath or reservoir. After the enormous irrelevant degrees of freedom of the thermal bath 
are integrated out from the von Nuemann equation of the density matrix of the extended system including the 
environment, the master equation for the reduced density matrix of the relevant system is resulted in some reasonable 
approximations; for details see, for example, Ref. @, 13)- F° r a simple system such as a two-level atom interacting 
with the background thermal radiation field, as the parameters of the master equation are independent of time, 
the solutions have been obtained and studied very well [USE!- However, for the non-autonomous case where the 
parameters of the master equation are dependent on time, Ellinas et al. have studied the adiabatic evolution and the 
corresponding Berry phase in optical resonance |9j , and the problem for the general solutions still remains open. For 
more complicated systems, the solution to the master equation of the reduced density matrix is, in general, difficult 
to obtain. To solve the problem, the usual way is to convert the master equation to a set of differential equations 
for some quantum statistical moments or expansion coefficients in terms of some bases truncated at a reasonable 
In the meanwhile, some other useful approximate methods have also been prop osed, for instance, the 
short-time expansion [Tfj| , the small lose rate expansion , the stochastic unraveling |12|. and finally numerical 
calculations 0] , etc. Some exact methods have been also explored 0, . An elegant method was proposed by 
Briegel and Englert to treat the quantum optical master equations by using the damping bases [l6| , but the problem 
was restricted to the autonomous case where the parameters of the master equation are independent of time. 

In this paper we shall present a novel algebraic dynamic method to solve the master equations which in many cases 
are found to have some dynamical algebraic structures. The common feature of the quantum master equations is 
the existence of the sandwich terms of the Liouville operators where the reduced density matrix of the system is in 
between some quantum excitation and de-excitation operators, which come from the elimination of the environment 
degrees of freedom. The sandwich structure of the quantum master equation also appears in thermal field theory 
where the so-called thermal Lie algebra has been proposed to treat the problem [ItI Il8j | , but the parameters of the 
master equation are still time- independent. 

Our new algebraic method is just a generalization of the algebraic dynamical method [l9| from quantum mechanical 
systems to quantum statistical systems with time-dependent parameters. It is designed to treat the sandwich problem 
and to explore the dynamical algebraic structures of the master equations with time-dependent parameters built in. 
To this end, the right and left representations |20j| as well as the adjoint representations of dynamical algebras 
are developed and used. This new method has been used successfully to solve the von Nuemann equation for the 
quantum statistical characteristic function of the two-level Janes-Cummings model |2l| and the master equation for 
the sympathetic cooling of the Bose-Einstein condensate system in the mean field approximation 22]. In this paper, we 
shall apply this method to solve the master equation of the dissipative two-level atom system in the non-autonomous 
case, and to study the decoherence of the multiple two level atom system coupled to the radiation thermal bath, 
which have not been studied before for the non-autonomous case. 
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The paper is organized as follows. In Sect. II, the model Hamiltonian of the system is given and the master 
equation for the reduced density matrix of the atom is described. In Sect. Ill, the dynamical u(2) algebraic structure 
of the Liouville or rate operator of the master equation is found by introducing the new composite algebras which 
are constructed from the right and left representations of the relevant algebras, and the dynamical symmetry of the 
system is established and the integrability of the master equation is thus proved by using the algebraic dynamical 
theorem 0|. Sect. IV is devoted to the analytical solutions of the master equation for the non-autonomous case 
where the parameters of the rate operator (or Hamiltonian) are dependent on time, and the approach to the steady 
solution asymptotically is thus proved. In Sect. V, the dissipation and decoherence of the multiple two-level atom 
systems are investigated for the non-autonomous case. Discussions and conclusions are given in the final section. 



II. DISSIPATIVE TWO-LEVEL ATOM SYSTEM IN NON-AUTONOMOUS CASE 



Consider the two-level atom interacting with a thermal quantum radiation field. With the dipole interaction and 
in the rotating wave approximation, the total system is described by the following model Hamiltonian, 

H = hiujocrz + Hbath + h{a+A + a-k 1 ), (1) 

where Hbath — X)fe J^kbtbk describing the background quantum radiation field, A = ^ fe gkbk is an operator used to 
describe the coupling between the atom and the radiation field with coupling constants g^. Here ojq and u>k are the 
transition frequency between two levels of the atom and the mode frequencies of the radiation field, respectively. er+ , 
er_ , and er 2 are dimensionless atomic operators obeying the usual Pauli matrix commutation relations |23| . Using the 
standard technique from quantum optics, one obtains the master equation for the reduced density matrix of the atom 

1 

= ~u [a z ,p{t)) - ~ ("o + 1) (o+cr_p(i) + p(t)<T + cr- - 2a-p(t)a + ) 

1 - 

-- n (a-a+p(t) + p{t)o-o + - 2a + p(t)a^) (2) 

where no is the mean number of photons in the environment and 7 denotes the damping rate. They read 

n = [exp(hu) /k B T) - l] -1 , 

7 = 2w^2gl5(w -u k ) (3) 

k 

Here the term which gives rise to a small Lamb frequency shift Aw has been neglected. Equation (J2J describes an 
atom interacting with a thermal field at the temperature T. If T = 0, then no = 0. 

In the autonomous case, the system has been studied very well. However, it is also interesting to control the 
system through changing the temperature of the thermal bath, the atomic energy level and the coupling constant. 
In this case, the parameters of the rate operator 7, no, and loq are time-dependent, and the system becomes non- 
autonomous. More basically, even if the total Hamiltonian of the composite system-the system to be investigated plus 
the environment, is autonomous, the master equation of the reduced density matrix of the investigated system still 
becomes non-autonomous under the non-Markovian dynamics |24j |. Therefore, the quantum master equation of the 
reduced density matrix, in general and more rigorously, should be non-autonomous in the sense that its parameters 
should be time-dependent. To our knowledge, the problem of the non-autonomous systems has not been solved till 
now. Therefore, in the present paper it is our main goal to solve the problem of the non-autonomous case for the 
two-level system. In the next section we first study the algebraic structure of the master equation© and explore its 
dynamical symmetry. 



III. ALGEBRAIC STRUCTURE OF THE MASTER EQUATION 



A. Right and left algebras in the von Neumann space 



Following the idea of Ref. |22|, the right and left algebras are introduced. First, one notices that the density matrix 
p is a super vector in the von Neumann space pflj . 



p = J2P* S '\s)(A (4) 

s,s' 
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where \s) denotes the Fermion state. The o + , er_, and a z can operate on the ket state |s) to the right and on the bra 
state (s| to the left, which form the right and left representations of the usual su(2) algebra as follows [2(| 

su(2) R :{a r z ,a r + ,aL}, su{2) L : {a[, a l + ,a l _}. (5) 

They obey the commutation relations respectively as follows, 

[a r z ,a r ± ] = ±2<4, [o r + ,o r _]=*l, 

[a l z ,a l ± ]= T 2a l ± , [a l + , a l _] = -u\. (6) 

It is evident that su(2)r is isomorphic to the su(2), while su(2)l is anti-isomorphic to the su(2). This is because that 
su(2)r operates, as usual, towards the right on |s). On the other hand, the su(2)l operates towards the left on (s\. 
Since su(2)r and su(2)l operate on different spaces(the ket and the bra spaces), they commute each other, i.e. 

[su(2) L ,su(2) R }=0. (7) 



B. Composite algebra and algebraic structure of the master equation 



After having introduced the left and right algebras in Eqs. (jSJ), one can construct the composite su(2) and u(l) 
algebras as follows, 

su(2) : {j Q = ^l±Aj + = a r + a l _,j^=aLa l + }. 

«(1) : Uo = (8) 
According to Eqs. © it is easy to check the following commutation relations 

[Jo, J±] = ±2J±, [j + ,J_]=J , 

[#o, 4] = 0, [U Jo] = 0. (9) 
The action of the composite su(2) and u(l) algebras on the bases of von-Neumann space is 

Jo\s)W\ - S -^\s){s'\, 

J+\s)(s'\ = 6 s+lfi 5 s , +lfi \s + 2)( s ' + 21 

J-\ S )(s'\ = ^-i,o^_in|s-2)( S , -2| ) 



where s and s' are equal to -1 or +1. 
Noticing the following identities 



Uo\s)(s'\ = S -^\s)(s'\. (10) 



cr + cr - = — ^ — ' a ~ a+ = — 2 — ' ^ ' 



and the composite algebra Eqs. (JSJ) introduced above, the master equation(0l can be rewritten as 

dp(t) 



dt 



T(t)p(tl (12) 



where the rate operator f reads 



f = -w (t)Uo +7(*)no(t)4+7(*)[«o(*) + 1]-^- - \l{t)Jo ~ ^W[2SoW + 1] (13) 

The master equation i|12|) is now put into a form similar to the time-dependent Schrodinger equation except that the 
imaginary number "z" is missing from the left hand side and the rate operator T is non-Hcrmitian, indicating the 
dissipative behavior of the system due to the energy exchange with the thermal bath. In Eq. 1)12(1. the rate operator 
f plays the role of the Hamiltonian and the reduced density matrix plays the role of the wave function. Since the 
rate operator T is a linear function of the su(2) ©m(1) generators, the master equation (|12fl possesses the su(2) ©w(l) 
dynamical symmetry; it is integrable and can be solved analytically according to the algebraic dynamics [19| . 
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IV. EXACT SOLUTION TO THE MASTER EQUATION IN NON-AUTONOMOUS CASE 

A. Eigensolutions of the rate operators T 

In order to better understand the time-dependent solution of the master equation, its decay behavior and its 
approach to the equilibrium state, we first consider the steady eigenvalue problem of the rate operator T whose 
eigensolution itself is interesting and peculiar, and contains the steady equilibrium state. The eigen equation reads 

tp(s,s')=(3(s,s')p(s,s'), (14) 

where [3{s, s') is the eigenvalue of the rate operator T and (s, s') label the eigenstates in the von Neumann space. This 
eigenvalue equation is time-independent and can be solved by introducing the following similarity transformation 

p(s,s') = Up(s,s'), (15) 

where 

U = e a + S +e a - J -, If- 1 = e -°>-j- e -<*+j+, (16) 

Here a± are the parameters specifying the similarity transformation. After some calculations, one has the transformed 
eigenvalue equation as follows, 

tp(s,s') = f3(s,s')p(s,s'), 

T = U- 1 fU 1 (17) 

where the transformed rate operator T is diagonalized and becomes a linear combination of the commuting invariant 
operators Jo and Uq which dictate the dynamical symmetry of the system, 

f = -iu U ~ ^7[2(n + l)a+ + 1] J - \l{2n + 1) (18) 
if the following diagonalization conditions are fulfilled 

— («o + 1)"+ - aq- + no = 0, 

(n + l)(l + 2a + a_)+a_ = 0. (19) 

The eigenvectors of f are the common solutions of Jo and Uo, just the form of s)(s' with s and s' = ±1. Eqs. i|19|) 
have two sets of solutions, which yield two sets of eigenvalues j3{s, s') for the rate operator 

r \ -, no + 1 „, s - s' 7,„_ ^..s + s' . 
(«) a+ = -l, a ~ = 2 fi a + V ^( s ' s ) = _ia; o— y~ + -^{ 2n + !)(— ^ 1), 

(b) « + -^T: /9 ( fl y) = - iwo i^ + 2(2fi + l)(-i±^-l). (20) 

no + 1 2no + 1 2 2 2 

At first glance, it is surprising that two similarity transformations exist to diagonalize the same rate operator and 
to yield two sets of eigenvalues. This is in contrast to the diagonalization of a Hamiltonian where the unitary 
transformation to diagonalize the Hamiltonian is usually unique and the set of the eigensolutions is also unique. The 
peculiar results stem from the special structure of the rate operator (|13fl : (1) it contains a part (from the second to 
the forth terms) which is a vector in the linear space spanned by the su(2) generators( J+, J_, Jo), and allows two 
transformations to rotate this part of vector along the Jo and — Jo directions; (2) the last term of the rate operator 
is a constant term which is a scalar in the su(2) space and make above two diagonalizing transformation asymmetric. 
The above two features result in two similarity transformations p^ . However, as will be seen soon, after returning 
to the physical frame by the inverse transformations, the two sets of eigensolutions coincide. This means that the 
physical results are objective, independent of the similarity transformations used. 

It is interesting to note that both solution (a) and (b) contain the zero-mode steady solution and the nonzero- mode 
decaying solutions( with negative eigenvalues ), which guarantee that any time-dependent solution of the reduced 
density matrix asymptotically approaches the steady solution as shown below. 

Performing an inverse transformation, the eigensolutions of the rate operator is obtained readily 

p{s,s') = Up(s,s') = (l + a + J + )(l + a_J_)|s)(s'|. (21) 
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Explicitly, both (a) and (b) solutions lead to the same physical eigensolutions. 

ft = /3(-l,-l) = 0, pi=p (_i,_i) = J£±L_i)<_i| + | + 1)<+1|, 

2no + 1 2n + 1 

02 = /3(+l,+l) = -7(2no + l), P2 = +1) = | - 1)(-1| - | + 1)(+1|, 

A = /?(+!, -1) = - J(2n + 1) - two, P3 = p(+l, -1) = I + 

/3 4 = /?(-l,+l) = -2(2no + l) + iw 0) p 4 =p(-l,+l) = | - 1>(+1|. (22) 

where the first line of Eqs. I)22H is the zero-mode solution corresponding to the steady state. 

Another feature of the rate operator is its non-Hermiticity, i.e., ft ^ f , which is evident from jt = J_, jl = J + , 
Jq = Jo, and Uq = Uq. Since f is non-Hermitian, the eigenvectors of f and ft constitute a bi-orthogonal basis pfij . 

Using a similarity transformation J7 = e ~ a + J - e ~ a - J + an(1 under the conditions of Eq. I|19|) . diagonalization of the 
operator ft can be obtained as follows 

ft = iuj Uo - ^7{2(n + l)a+ + 1]Jq - ^7(2n + 1) (23) 



In a similar way, we get the eigensolutions of ft . 

PI = /3*(-l,-l) = 0, p 1=i 5(-l,-l) = |-l)(-l| + | + l)(+l|, 

/?* = ^(+i J+ i) = - 7 (2fio + i), p 2 = p(+i,+i) = -^|-i)(-i|-|^-|--i);h-i|. 

znp + 1 zno + 1 



(it = f3*(+l, -1) = -2(2fio + 1) + two, p 3 = P(+1, -1) = I + 

/3 4 * = J 9*(-i,+i) = -l(2n + l)-iw , p 4 = p(-l,+l)=|-l)(+l|. (24) 

where (3* and /5j are the eigenvalues and eigenvectors of the operator ft. It is can be checked that Eqs. 1221) and Eqs. 
(|24|l are bi-orthogonal. 

B. Time-dependent solutions of the master equation in Non-autonomous case 

Since Eq. Q13f) is a linear combination of the generators of the composite algebra, the master equation Eq. I|12|) 
possesses the sm(2)©u(1) dynamical symmetry and is thus integrable, solvable analytically even in the non-autonomous 
case [l9|. 

The master equation in the non-autonomous case can be solved by algebraic dynamical method via the following 
time-dependent gauge transformation which is a generalization of time-independent similarity transformations in the 
autonomous case to the non-autonomous case ( the terminology of gauge transformation is due to the fact that it 
induces a gauge term in the rate operator similar to the gauge field theory), 

U g = e a ^ t)J +e a - {t)J -. (25) 
After the gauge transformation and under the following gauge conditions 

da+(t) . 2 . . , . . , . no(t) 



-7(*)[no(*) + l]<(t) - l(t)a+(t) + 7(i)no(t) = -7(*)fa>(*) + !][«+(*) + !][«+(*) 



^=7(()[n„(()|l][l|2Q + (()L(()]+ 7 (() ft _(i), (26) 
the gauged rate operator is diagonalized and the gauged master equation becomes simple and integrable, 

dm mm, 



df 

T(t) = U-'mUg - U; 1 ^- = -iuo(t)U - ^7(*){2[n (t) + l]a+(t) + 1}J - \l{t)[2n {t) + 1]}. (27) 
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which is a linear function of the complete set of the commuting operators (invariant operators)[/o and Jo, and clearly 
shows a u(2) dynamical symmetry. The solution of Eqs. (|27[) reads 



For the initial conditions 



p(t)=e/o r W%( ), (28) 

o+(0) = a_(Q) = 0, 
orp(0) = p(0)=5^p„/|a)<a'| (29) 

ss' 

or , we finally obtain the solution, 

p(t) = e a +( t )^+ e a -( t )^- e ~ i ^° fg ^o(r)dT-J f o t 7(r)((n (T) + l)Q + (r) + i)ciT- J* 2^1 (2fl (i-)+l)dT /^qx 

Once the reduced density matrix of the non-autonomous system is obtained, the averages of the physical observables 
<r z , cr+, and cr+ can be calculated. For the system initially in a pure state, we have p(0) = |yLt| 2 |l)(l| + |^| 2 | — 1 + 
p,v*\l)(— 1| + p?v\ - 1)(1|, where \p\ 2 + \v\ 2 = 1. Then one obtains 

Pit) = [/i,iWIm| 2 (1 + «+(*)a-(t)) + /_ 1) -i(t)MV(*)]|l><l| + [/ M (t)|u| 2 a_ (t) + /-i,-i(t)M 2 ]| - 1)(-1| 

+/x,-i(t)/i^|i>(-i| + /. w (iKi/| - (3i) 

where 

/„ s ,(f) = e - l£ T £l X t " ( '( T ) dT - £ T £: J t 7(^){[fio(r) + l] Q+ (r) + i}rfT-i 7 (r) [2n (r) + l]<ir ^ 

From Eq. (|31fl. we get 

(<r»> = /i,i(t)|M| 2 [l + a+(t)a_(t)-a_(t)] + /_i,_i(t)|z/| 2 [Q! + (t)-l], 
(o+> = f-i,i{t)p*v 

(<r_) = /i,_i(t)/«>*. (33) 
For the autonomous case, no and 7 are independent of time, a+ and a_ have the following analytical solutions, 

1 _ e -7(2fio + l)i 



«+(*) 



"0 + 1 _|_ e -7(2n +l)t : 



(no + l)n [^Ji +e-^( 2 " 0+1 ) t ][l - e -r(a«o+i}*] 

= (2n + l)2 e -7(2n +l) f ' ^ 



Then 



(2n + l)e'^( 2 "" +1 )* 
(no + 1) + n e-T(2*io+i)t' 

(n + 1) + n e- 7 ( 2fio+1 ' t 



2n + 1 
/l -l(*) = e _<Wo * _ 3K 2 ™ 0+1 ) t , 

/_i,i(t) = e ^o*-i( 2 »o+i)*. (35) 
Inserting Eqs. (|34|l and Eqs. (|35|l into Eqs. (|33|l . we recover the well-known results for the autonomous system. 



C. Approaching the steady solution 

Like the autonomous case, the time-dependent solutions of the master equation in the non-autonomous case also 
asymptotically approaches the steady solution satisfying 



(36) 
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which has the same solution as the zero-mode eigensolutions of T , namely the steady solutions with the parameters 
a + and a_ obeying 



da+(t) = Q 

dt 
da-{t) 



dt 



0. (37) 



whose solutions are obviously the same as Eqs. i|19|) and should have the two same sets of solutions (a) and (b) as given 
in Eqs. 1)20(1 . For the autonomous case, the two similarity transformations as tools to diagonalize the rate operator 
are on the equal footing and generate the same physical solution as proved above; while for the non-autonomous 
case, the properties of the two sets of solutions of a+ and a_ should be examined on the background of their time 
evolution. It is found from the Eq. (26) that if a + = — 1 — e,then da ^^ < and a+ will go away further from — 1 

towards negative direction; while as a + = -1 + e, then > q an d a+ will go away further from —1 towards 

positive direction. Thus the steady solution a + = — 1 is unstable. Since the solution of cx-(t) depends on a + (t), it is 
also unstable. Therefore the solution (a) is unstable and can not be reached from the initial condition a+ = a_ = 0. 
Instead, the solution (b) is partly stable in the following sense. Because the first equation of Eq. (26) reads 

(! " - </;) = -7(*)[no(t) + l]Mt) + !][«+(<) - J^-}- (38) 



dt ,WL uv ' JL ^ w " ^ v ' n (t) + V 



It's easy to see that 



da+(t)/dt > 0, if < da+(t) < 0^ 



da+{t)/dt < 0, if a+(t) > and a+(t) 



< -1 



da+(t)/dt > 0(< 0) if da + (t)/dt > (if a + (t) > j^rfj and a+(t) < -1). With the initial condition a + = 0, we see 

»o(oo) _ ii 
fi (oo) + l no + 1 



that a + approaches the value -"f^f^ = -"^ asymptotically from zero. However, a_ (t) can not reach its steady value 



■>„ 0+ j- To study the asymptotic behavior of a_, we define y(t) = a_(f) x exp{— J Q 7(r)(n (r) + l)(a+(r) + l)dr} = 

a_(i)exp J^p(r)dT. The time differential of y (t) is given by frexp f Q p(r)dT where b = a~(t) + ct-(t)p(t). Since 
b — ► 7(^0 + 1) which is bounded and p(t) is negative for large t, the differential tends to zero and, hence, y (t) is 
towards a constant. This implies that a~(t) diverges asymptotically. So the unique solution of Eq. (26) has the 
following asymptotic properties: 

a_(t) X e -/ S«(no(r)+l)( a+ (r)+l)dr (3g) 

Using the above asymptotic relations, one obtains the asymptotic results of the time-dependent solutions as follows 

P+-{t) = e"^o" o(r)dT -/o^ (T)( " o(T)+ 5) d V +(t)J +e Q -W J i + l)(-l| 

= g-i/o w o(T)dT-/ o ' 7 (r)(fio(r)+i)dT| + y ^ 

*>-+(«) |t-oo = e +i /o Uo(T)iT -/o 7(T)(fio(T)+ ' )liT e a +( t )-'+e a -( i !''-| 

|t-oo = e "/o^ r)( " o(T)+1)(Q + (T)+1)dT e Q +( t ) J + e Q -( t ) J i + l)(+l| 

= e - r o , 7 (r)(ft (r)+l)(a + (r)+l) d r [| + + + q _^ )(| _ + a+(f) | + 1)(+1 |Vj 
► Const. X (I - 1> < — 1 1 + _ — I + 1)(+1|) = Const. X p st eady, 

n + 1 

p (£) = e -/ *7(T){-(fto(r) + l))a + (r)+fi (T)}dr eQ+ ( t )J +eQ „( t )J_| _ 

= e -/ t 7(r)[-(ao(r)+l)a + (r)-H i o(r)]«»r [ | _ X) + a+(t) | + _^ c x (4Q) 
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In the above derivation, we have used the following asymptotic relations: 

e ~ J 7(T)[-(«o(-r) + l)(a + (r)+fio(T)]dr ^ ^ 

which can be proved by Eq. (26). 

The above results indicate that the time-dependent solutions of the master equation in the non-autonomous case also 
asymptotically approach the steady solution irrespective of their initial conditions, and that the divergent behavior 
of a>-(t) and eq.(37) play an important role in the process of approaching equilibrium state. 



V. DECOHERENCE OF MULTIPLE ATOM SYSTEMS 



It is straightforward to generalize the above model to N two-level atoms. The problem of the N two-level atoms 
system is related to the quantum register and the entanglement state in quantum computation |26j ] . Because of the 
inevitable coupling of the atoms to the external environment, the entanglement state will lose the coherence among 
different atomic states and some information carried by the multiple atoms will be lost. Palma et al. has studied the 
impact of decoherence on the efficiency of the Shor quantum algorithm and the decohcrcncc of quantum register at 
the two qubit level [2?|]. However, the solution to the problem in the non-autonomous case is still lacking. The N 
two-level atoms coupled to the quantum radiation environment can be described by the following model Hamiltonian, 

H = H s + H env + Hi, 
1 N 
fc=i 

N oo 

fc=U=i 

Taking the same procedure and the same approximations as those for the one atom case , the master equation for the 
N two-level atoms can be obtained 

N 

f>N = ^2 = r PJV 

k=l 
N 

r = £r fc , 

k=l 

T k = -twot/f + 7 «oi++ l(n + l)Jl~\ Jo -|(2no + l). (43) 

Here we have assumed that the N two level atoms are identical and coupled to the same environment so that 
the decay rate 7 for different atoms and the mean number of environment photons no are the same (if they are 
different, a superscript "k" should be put on each pair of parameters, namely 7 fc and ). We have also set 
fli = H2 = ■ ■ ■ = £In = ^o- The solution to Eqs. I|43l) reads 

N 

PN(t) = l[p k (t). (44) 

k=l 

According to Eq. (|30(l . the density matrix pk(t) of the k-th qubit is 

Pk{t ) = e<W J ie»-W J -eJo tk ^ Pk (0), 

T k (t) = -iu; Q (t)U k - i 7 (*){2[no(t) + + l}^o 

-I 7 (t)[2n (t) + 1]. (45) 

where T k (t) = -iu (t)U$ - ±j(t){2[n (t) + l)a k + (t) + l}J fc - | 7 (i)[2n (t) + 1]. As in the last section, the initial 
density matrix pfc(0) of each qubit can be expanded in terms of the superbases, pk(0) = Y^u s') c s,s' \s) kk (s'\. Here 
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a^_(t) and a'l(i) obey the Eqs. because we have assumed that the parameters of the rate operators for different 

atoms are the same functions of time. If ^ k and n§ are different for different atoms, aq_(f) and a_(t) are also different 
for different atoms. However they obey the same form of Eqs. (|26[1 but with different parameters 7 fc and Uq. To 
illuminate the model concretely, we also consider the two qubit system with the initial state in the pure state, namely 
|^(0)) =a\ + -)+/3\ -+), and 

p 2 (0) = |a| 2 | + l) n (+l| <g> | - I) 22 (-l| + |/3| 2 | - 1) U {-1| ® I + 1) 22 (+1| +af3*\ + l) n (-l| <8> | - I) 22 (+l| 

+a*/3\ -l) n (+l|® | + I) 22 (-f| (46) 

The time-dependent solution of the density matrix is now 

P 2(t) = |a| 2 /i,i(*)[(i + «+(*)«-(*))! + i) n (+i| + «-(*)! - i) u (-i|] ® (*) [I - i) 22 (-i| + «+(*)! + i) 22 (+i|] 

+ |/3| 2 /_i,_i (i) [| - l) n (-l| + a+(t)| + I) n (+l|] ® A,i(t)[(l + a+(i)a_(t))| + I) 22 (+l| + a_(t)| - 1> 22 <— 1|] 
W/wl + l) n (-l| ® - 1) 22 (+I| + - l) n (+l| ® + 1} 22 {-1| (47) 

where f SlS i (t) are the same as given in Eqs. <|32|) . For the autonomous case, 7(f), no (f) & n d wo(f) are independent 
of time, and the solutions of cx+(t) and a_(f) are determined from Eqs. i|26|) . which are the same as Eqs. (|34|) . Eq. 
(14711 is now reduced to the result for the autonomous case, it reads 



_ e -2 7 (2, 0+1)t4 Mno + l) p x ^ ^ + p} p2 (48) 

The above solution indicates that during the time evolution, the density matrix of the two qubit entanglement state 
will approach the steady density matrix( the last term in Eq. I|48|l )and lose its coherence. The characteristic time of 
the decoherence is T decoh = 7 ( 2fi 1 0+1) ■ 



VI. SUMMARY 



In this paper, Based on the quantum master equation, we have investigated the dissipative and decoherence behav- 
iors of the two-level atom system coupled to the environment of thermal quantum radiation fields. The dynamical u(2) 
algebraic structure of the quantum master equation of the two-level dissipative system in the non-autonomous case is 
found by virtue of left and right algebras. By the algebraic dynamical method and proper gauge transformations, the 
analytical solutions to the non-autonomous master equation are obtained and the long time behavior of the system has 
been examined. Finally we extended the model to the multiple two-level dissipative atom system and its decoherence 
is studied in terms of the density matrices for the non-autonomous case, which are given analytically and related 
to quantum register and quantum computation. Since the master equations of a wide class of dissipative quantum 
systems possess some dynamical algebraic structures, the present method used by us may serve as a useful tool in 
quantum statistical physics to treat the dissipative and decoherence problems. In addition, the results obtained in 
this paper may be practically useful for the analysis of the decoherence of the multiple two-level atom systems and 
quantum register. 
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